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The present work provides the results of a short exploratory study on the performance 
of Phenolic Impregnated Carbon Ablator, or PICA, at high heat flux and pressure in an 
arcjet facility at NASA Ames Research Center. The primary objective of the study was to 
explore the thermal response of PICA at cold-wall heat fluxes well in excess of 1500 
W/cm. Based on the results of a series of flow simulations, multiple PICA samples were 
tested at an estimated cold wall heat flux and stagnation pressure of 1800 W/cm/2 and 130 
kPa, respectively. All samples survived the test, and no failure was observed either during 
or after the exposure. The results indicate that PICA has a potential to perform well at 
environments with significantly higher heat flux and pressure than it has currently been 


flown. 
Nomenclature 
CFD = Computational Fluid Dynamics Dp = Base diameter 
CT = Computed Tomography Abulk = Bulk total enthalpy 
DPLR = Data Parallel Line Relaxation code Aci = Centerline total enthalpy 
FIAT = Fully Implicit Ablation Thermal K = Constant of proportionality ina 
response code simplified Fay-Riddell correlation 
IHF = Interaction Heating Facility Dstag = Stagnation point pressure 
IM = Inspiration Mars stag = Stagnation point heat flux 
IR = Infrared Rp = Base radius 
MSL = Mars Science Laboratory Re = Corner (or shoulder) radius 
MSR = Mars Sample Return Regt = Effective radius 
PICA = Phenolic Impregnated Carbon Rn = Nose radius 
Ablator 
TPS = Thermal Protection System 


I. Introduction 


hermal protection systems (TPS) for capsules that execute atmospheric entries at hyperbolic velocities (in excess 
of 11 km/s at Earth) are based on ablative materials. One class of ablative materials is low density, Phenolic 
Impregnated Carbon Ablator, or PICA, which is of particular interest, because its low density implies a relatively light 
TPS or heatshield, and hence, a lower entry mass of the capsule. During entry heating, the phenolic phase pyrolyzes 
and the pyrolysis gas blows the hot boundary layer away from the wall. A charred region composed of carbon is left 
behind after pyrolysis is complete. Any further heating causes the carbon to either oxidize, erode or sublimate, with 
consequent surface recession. 
PICA was originally developed for the NASA Stardust mission, which successfully returned samples of comet 
material back to Earth in 2006.' During reentry of the Stardust capsule at roughly 13 km/s the PICA heatshield (built 


' Senior Research Scientist, Thermal Protection Materials Branch, Mail Stop 234-1. Senior Member AIAA. 
* Senior Staff Scientist, Entry Systems and Technology Division, Mail Stop 229-1. Associate Fellow AIAA. 
> Aerospace Engineer, Thermal Protection Materials Branch, Mail Stop 234-1. Senior Member AIAA. 

+ Senior Research Scientist, Thermal Protection Materials Branch, Mail Stop 234-1. 


1 
American Institute of Aeronautics and Astronautics 


in a monolithic construction) was exposed to an estimated peak heat flux of 1000 W/cm7, and a stagnation pressure of 
28 kPa (0.28 atm).” A slightly different formulation of PICA was used as the heatshield (in a tiled construction) as 
shown in Figure 1, for the Mars Science Laboratory (MSL) spacecraft, which carried the Curiosity rover to Mars.° 
MSL did not experience the same level of heating as Stardust. PICA is currently the baseline heatshield material for 
the OSIRIS—REx sample return mission* scheduled to launch in 2016. Furthermore, SpaceX uses a variant called 
PICA-X on its Dragon spacecraft that ferries cargo to and from the International Space Station.° 

Since the successful demonstration of PICA’s flight performance up to 1000 W/cm? and 0.28 atm of pressure, 
several subsequent successful tests in the arcjet facility at NASA Ames 
Research Center (ARC) have shown PICA material to be capable of 
handling heat fluxes up to 1600 W/cm7, at surface pressures up to 130 kPa.° 
If indeed PICA would exhibit no failure at heat fluxes and pressures higher 
than previously measured, it could become viable as a heatshield material 
for other missions, e.g., the Inspiration Mars (IM) mission concept of 
Dennis Tito as well as Comet Surface Sample Return (CSSR) missions. 
The Inspiration Mars mission concept explored the feasibility of sending 
two people on a fast fly-by of Mars in 2018 with a propulsion-free return 
to Earth at roughly 14 km/s.’ Results of the mission study showed that the 
entry capsule would experience a peak heat flux of the order of 2000 
W/cm”. The questions then are: (1) can standard PICA be tested at higher 
than 1500 W/cm? heatflux in the arcjets and (2) in the event of failure of 
the standard material, could it be densified — at the surface or through the Figure 1. MSL PICA heatshield 
thickness — to extract ablative performance without material failure? 

Toward this end, tests were performed in the arc-heated flow through a 15 cm (6 inches) nozzle of the 60 MW 
Interaction Heating Facility (HF) at NASA ARC. The primary objective of these tests was to explore the performance 
of standard PICA and its denser variants at high heat fluxes (in excess of 1500 W/cm?) and high pressures (in excess 
of 1 atm). The emphasis of the present paper is solely on the performance of standard PICA. The sections that follow 
provide details of the test design and the key findings from the test. 





II. Arc-Jet Model Design and Pre-test Aerothermal Analysis 


The overall geometry of an arcjet stagnation test article consists of two major parts — a leading portion made of the 
TPS material to be tested, and a model holder made of graphite; the model holder is attached to the water-cooled sting 
of the arc jet facility. In order to optimize the design for test sample and achieve the desired test environments pre-test 
aerothermal analysis was performed. The test sample (TPS material) diameter was set to 5 cm (2 inches) to achieve 
heat flux values > 1500 W /cm’. While we could have used smaller samples to achieve even higher heatflux, the 
sidewall heating would have been significant and for the present tests we decided to go with 5 cm diameter. The holder 
design and material selection was also critical to avoid any undesired failure events (e.g. sample falling off the holder) 
either during or immediately after the plasma exposure. Figure 2 shows the various axisymmetric geometries evaluated 
computationally as part of the aerothermal design phase. The three categories of geometrical considerations in the 
aerothermal design are shown in Table 1. 
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Figure 2. Candidate axisymmetric geometries evaluated in the aerothermal design phase. The overall geometry consists 
of two pieces — a leading portion made of the material to be tested, and a model holder made of graphite. 
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Table 1. Geometrical considerations in aerothermal design of the test article 


Primary Design Variables | Comments 


TPS material Radius of curvature of the | Curvature is required to keep the heat flux high 

exposed face at the stagnation point; a flat-faced article will 
experience lowest heating at the stagnation 
point 

Shoulder radius The radius of the shoulder has an influence on 

the stagnation point heat flux (work of Zoby 

and Sulllivan®); the smaller the shoulder 
radius, the lower the stagnation point heating. 

Angle of the side wall This variable is driven by sidewall heating and 

consequent influence on any embedded 

thermocouple; a sidewall angle of 0° (cylinder) 
yields a low heat flux distribution along the 
wetted length. 

Length of the test article This variable is driven by the requirement that 

there be enough virgin material left to be able 

to clearly discern char and density change over 
the depth (thickness) of the material. 

Graphite holder This variable is fixed by the test facility. 

Ramp angle The taper angle of the graphite holder is 
determined by the level of axial flow 
separation at the foot of the ramp; although 
flow separation provides some heating relief 
over the length of the separated flow region, 
the holder will experience locally enhanced 
heating at the point of impingement of the 
shear layer. 

Test Material-Holder interface | Step height A backward facing step is necessary to force 
axial flow separation; the step height is 
determined by the size of the separation 
bubble. Furthermore, a small gap between the 
test material and graphite holder is necessary 
to account for disparity in thermal expansion of 
the two parts. The size of the gap is driven by 
thermal-structural considerations. 








As mentioned earlier, 15 cm (6 inches) exit diameter nozzle of the 60 MW IHF was chosen to achieve highest 
possible heating environment on the test article with that nozzle. For heater settings of 6000 A of current and a mass 
flow rate of 0.85 kg/s, nozzle computations were performed using the analysis procedure developed by Prabhu et al.’ 
This procedure uses the NASA flow solver, DPLR,'° as the main computational engine. DPLR solves the equations 
governing the flow of a gas mixture (N2, O2, NO, N, O, and Ar for the IHF) in thermochemical nonequilibrium. From 
the nozzle free jet solution, flow variables — mass density, velocity, temperatures (translational and vibrational), and 
mass fractions of the component species — were extracted at a point on the nozzle axis and 10 cm (4 inches) from the 
exit plane of the nozzle. These flow variables were used as uniform freestream conditions for the various candidate 
test article geometries. Computations for the test article were performed with two different wall boundary conditions: 
(a) wall temperature fixed at 400 K to mimic a calorimeter simulation, and (b) wall temperature allowed to float with 
imposition of radiative equilibrium (albeit without conduction through the thickness). 

Contours of Mach number in the pitch plane of each of the six candidate geometries are shown in Figure 3. The 
contour plots provide a qualitative picture of the flow. The plots clearly show the influence of the ramp angle of the 
graphite holder — the steeper the ramp angle (Baseline and Design #3), the larger the region of axially separated flow. 
A shallow ramp for the model holder ensures the flow stays attached on the test coupon, and increases sidewall heating. 
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The surface pressure and heat flux distributions (Figure 4) for each of the configurations give a quantitative picture 
of the influence of various geometric parameters of the test models. Results of these pre-test computations showed 
that it might be possible to achieve a cold-wall heat flux of 2000 W/cm? (at a pressure of 1.4 atm) at the stagnation 
point of the test article. This predicted heat flux is a 25% increase over the cold-wall value of 1600 W/cm? (at a 
pressure of 1.2 atm) reported by Beerman and Stackpoole.°® 
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Figure 3. Contours of Mach number in the pitch plane of candidate geometries. Axial flow separation is the 
largest for the steeper taper angles of the graphite holder. 
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Figure 4. Surface distributions of: (a) pressure, and (b) cold-wall heat flux for the 6 candidate geometries 
considered in the aerothermal design phase. Two of the configurations, designated as Baseline and Design #3, 
have significant axial flow separation at the foot of ramp and the high heating (= 400 W/cm2) over the graphite 
holder. 

The candidate geometry designated as Design #4 (see Figure 3) was selected as the final shape to be tested, and 
several models were constructed for tests in the IHF. The test coupon of 5 cm (2 inches) cylindrical diameter had a 
truncated spherical section with a radius of curvature of 5 cm (2 inches) and a shoulder radius of 0.64 cm (0.25 inch); 
an 1so-q body has its frontal radius of curvature the same as its base diameter, i.e., R, = Dp. The entire test coupon was 
6.4 cm (2.5 inches) long. The graphite holder had a taper angle of 15.3°, and at the test coupon-graphite holder interface 
there was a small rearward facing step of 1.3 mm (50 mil). The test coupon was attached to a graphite holder with 
help of 3 graphite pins and Graphi-Bond 551-RN (a high temperature graphite adhesive) at the interface, and the 
graphite holder, in turn, was connected to a water-cooled copper back plate to prevent the graphite holder from 
overheating. The overall length of the test articles ((ncluding the holder) was approximately 20 cm (7.8 inches). A 
dimensioned section of the entire test article and its interface to the facility sting are shown in Figure 5. The assembly 
of the various components of the test article is shown in Figure 6. 
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Figure 5: Model dimensions (cm) and pyrometer locations. 
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Figure 6: Components and assembly of test Article 


III. Pre-test Thermal Response Analysis 


The predicted aerothermal environments — surface pressure, heat transfer coefficient (heat flux divided by the 
difference between the edge and wall enthalpies), and recovery enthalpy — were then used as boundary conditions in 
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material thermal response computations code, FIAT,'! to estimate the time of exposure (to the arc-heated stream) 
required to recess the test material by 2.5 cm (1 inch) at the stagnation point. Surface thermochemical interactions and 
blowing effects were incorporated in the material-response code by use of ablation tables or dynamic chemistry, a 
surface energy balance with heat transfer coefficient, and a blowing reduction parameter. Further details on this 
coupling methodology may be found in the work of Milos and Chen. '” 


It should be noted that for most arcjet tests with an iso-q shaped sample, the surface environments predicted by 
DPLR are assumed constant over the time of exposure. However, if shape change or total recession are large, then 
the surface pressure and heat flux values and distributions may change significantly, in which case a second DPLR 
solution for the ablated shape at the recessed distance from the nozzle exit is also recommended. The initial and final 
environments then are interpolated in time for the material response computations. A different boundary condition, 
which allows radiation to the environment but no ablation, is used for the cool-down period after the test article is 
taken out of the arc-heated stream. While there is no imposed heat flux on the test article, the internal processes of 
heat conduction and pyrolysis continue during the cool-down period. 

The FIAT calculations for pre-test predictions for this test were performed using the nominal aerothermal 
environment, and also with a +10% scaling factor applied to the heat transfer coefficient. This factor is considered to 
be the minimum uncertainty of the arcjet environment. Figure 7 presents the predicted stagnation point recession 
histories for these three heating levels for the pre-test environment from Figure 4. At high heat flux, the recession rate 
of PICA is relatively steady. Based on these predictions, an exposure of duration of 20 seconds was proposed for the 
test series. 
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Figure 7. Predicted recession for pre-test environment with +10% variation in heat transfer coefficient. 


IV. Test Environments and Instrumentation 


The nominal arcjet conditions, based on the pre-test analysis and the historical data collected on the facility parameters, 
are listed in Table 2. The maximum possible heat flux on the sample was obtained by operating the arc heater at near 
maximum current and flow rates. The real time measurement of the centerline heat flux values was attempted by using 
51 mm iso-q calorimeters (identical to test samples) and 102 mm (4-inch) hemispherical calorimeter. The stagnation 
and sidewall surface temperature of the test samples were measured using optical pyrometers. Three different 
pyrometers were used for temperature measurements. The positions of pyrometers are shown in Figure 5. The first 
pyrometer was focused on the stagnation surface. The second pyrometer was focused 1.27 cm (0.5 inches) away from 
the stagnation surface to the depth where surface densification ends. The third pyrometer was placed 2.5 cm away 
from the stagnation surface, to get measurements at the sidewall. Surface temperature measurements using a fiber- 
coupled optical-emission spectrometer was also attempted during the test. In addition, an infra-red video was used to 
observe sidewall and stagnation surface temperature of the model. High definition video cameras were used to acquire 
videos during the arcjet exposure at various locations of the test samples. Pre- and post-test photographs of the test 
articles were also obtained. 

The measured stagnation pressure and heat flux can be used to infer the total enthalpy along the central streamline 
of the flow using a simplified version of the Fay-Riddell correlation.'* The simplified correlation for a cold wall is: 


 Psta 
Qstag = KE He (1) 
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where the constant K, which depends on the mixture composition, is taken from the work of Zoby,'* and Resr is the 
effective hemispherical radius of a non-hemispherical blunt body. The effective radius depends on the ratios R,/R, and 
R./Rp, where Rp, Rn, and R. are the base, nose, and shoulder (corner) radius, respectively, of the test coupon, and the 
value is determined from the work of Zoby and Sullivan.® In the present work, the effective radius of the 51 mm (2 
inches) 1so-q calorimeter is 45.4 mm; the effective radius of a 102 mm (dia) hemispherical calorimeter is the same as 
its geometric radius, i.e., 51 mm. 

The simplified correlation in Eq. | can be used to determine the consistency of heat measurements made by two 
different calorimeters at the same flow condition. Applying Eq. 1 to the hemispherical and iso-q calorimeters, and 
taking the ratio of the two expressions, one has: 

YUstag,iso—q = VRhemi (2) 
dstag,hemi JReff.iso—q 

Using the radii of the two calorimeters in Eq. 1, the ratio of heat fluxes is 1.06, i.e., the 51 mm iso-q calorimeter 
will measure 6% higher heat flux than the 102 mm hemispherical calorimeter at the same flow condition, a useful 
result that will be used later in post-test analysis. 


Table 2: Target Arc-jet environments for test 


Facility Nozzle Cold-Wall Stagnation Centerline Main Add = Argon Arc Distance 
Heat Flux Pressure Enthalpy Air Air Current from 
(2-1n 1sO-q) Nozzle Exit 
(cm) (W/cm?) (kPa) (MJ/kg) (g/s) (g/s) (g/s) (A) (cm) 
THF 15.24 2000 135 31.7 720 5D 52 6000 10.16 


V. Test Results 


A two day test campaign, consisting of 4 total runs (two each day), was successfully conducted in IHF. Calorimeters 
were included in the first run (day 1) and the third run (day 2). Three different MSL grade standard PICA samples 
were tested along with several surface-densified and high-density PICA samples. The results and post-test analyses 
focusing solely on standard PICA samples are described in the following subsections. 


A. Calorimetry and post-test CFD analysis 

The as-run heater settings and measured heat fluxes and pressures of the test entry (IHF 290) are listed in Table 3. 
The measured stagnation pressures from two different runs are in good agreement with each other. However, the 
measured heat fluxes do not appear to be consistent with each other. Melting of slugs was observed during the 
exposure, which could have been the cause of this inconsistency. 


Table 3: Heater settings and calorimetry data from test entry IHF 290 


Arc Main Add _ Argon Arc Bulk Calorimeter Cold- Stagnation Centerline 
Current Air Air column ~ enthalpy Wall Pressure Enthalpy 
pressure Heat Flux 
(A) (g/s)  (g/s)  (g/s) (kPa) (MJ/kg) (dia, geom) (W/cm) (kPa) (MJ/kg) 
6066 720 55 53 820 14.2 102 mm, 3259 130.9 50.5 
hemi 
51 mm, 2611 130.0 29) 
1SO-q 
6056 720 59 53 823 25 51 mm, 3507 131.7 47.3 
1SO-q 


The computational methodology used in the aerothermal design phase was applied to the new heater settings, and 
the resulting distributions are shown in Figures 8 and 9 for the hemispherical and iso-q calorimeters, respectively. 
Also shown in the figures are the corresponding measured data. 
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Figure 8. Surface distributions of: (a) pressure, and (b) cold-wall heat flux from DPLR simulations for the 102 mm (dia) 
hemispherical slug calorimeter used in IHF 290 (Table 3). Uncertainties of +5% and +15% have been assumed for 
pressure and heat flux measurements, respectively. 
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Figure 9. Surface distributions of: (a) pressure, and (b) cold-wall heat flux from DPLR simulations for the 51 mm (dia) 
iso-q slug calorimeter used in IHF 290 (Table 3). Uncertainties of +5% and +15% have been assumed for pressure and 
heat flux measurements, respectively. 


Although the agreement with pressure measurements is good for both calorimeters, the predicted heat fluxes are 
significantly below measurements. In order to verify the computational process, which has been anchored to several 
cases in the historical database of measurements in the IHF, an additional case was selected from a previous test, [HF 
216 with similar heater settings as the present test. The heater setting and calorimetry data from one run of this entry 
are shown in Table 4. The heater setting is identical to the one used in the aerothermal design phase. The predicted 
distributions of pressure and cold wall heat flux are shown in Figure 13. There is excellent agreement between 
measurement and prediction. 


Table 4: Additional heater setting and calorimetry data from test entry IHF 216 
Arc Main Add Argon Arccolumn Calorimeter Cold-Wall Stagnation Centerline 


Current Air Air pressure Heat Flux Pressure Enthalpy 
(A) (g/s) (g/s) (g/s) (kPa) (dia, geom) (W/cm7) (kPa) (MJ/kg) 
5985 740 55 54 868 102 mm, 1991, 152.1, 31.0, 
hemi 2090 132.8 32.5 
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Figure 10. Surface distributions of: (a) pressure, and (b) cold-wall heat flux from DPLR simulations for the 102 mm (dia) 
hemispherical slug calorimeter used in IHF 216 (Table 4). Uncertainties of +5% and +15% have been assumed for 
pressure and heat flux measurements, respectively. 


Comparing the results for the hemisphere (Figure 10) with those from preliminary design computations for the 
iso-q shape (Figure 4), the heat fluxes are consistent with the scaling in Eq. 2, which gave us confidence in the results 
of the computational method. The final test of the CFD results lies in the computation of the thermal response of the 
material to the predicted heat flux distribution. 


B. Performance of standard PICA samples 

Three MSL-grade PICA samples were tested in different runs. Figure 11 shows the top view and side view images 
of one of the PICA test samples prior to testing. Based on the calculations and recommendation from the FIAT thermal 
response model, the samples were exposed for approximately 20 seconds. Figure 12 shows one of the samples during 
plasma exposure. The sample was kept about 10.16 cm (4-inches) away from the nozzle exit. At approximately 20 
seconds the sample was removed from the plasma-stream. The top view and side view images of the sample after the 
exposure are shown in Figure 13. All the standard PICA samples performed very well during the arcjet plasma 
exposure. There was no evidence of spallation or any other damage during the test. Pre-test and post-test computed 
tomography (CT) scans were obtained for some of the samples. The images of a pristine (untested) PICA samples as 
well as arcjet tested samples are shown in Figure 14. It is evident from the CT scan images that ablation is very uniform 
and controlled in the PICA material. No unusual or concerning features or cracks were observed at the ablated surface 
or internally in the PICA material. 





Figure 11: Top view and side view of PICA samples prior to arcjet test. 
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Figure 12. Top view and side view of PICA samples during the arcjet test 








IHF 29 
Figure 13: Top view and side view of PICA samples after 20 second exposure to arcjet plasma. 





Figure 14: Pre-test and post-test CT scan images. 
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C. Surface temperature measurements 

Measurements from optical pyrometers for temperature measurements at the stagnation surface and sidewall were 
obtained for all three standard PICA samples. The data for the first PICA sample is shown in Figure 15. The stagnation 
surface turned out to be significantly hotter than the measurement limit of pyrometer and it saturated at 3000°C. FIAT 
simulations predict values > 3300°C at stagnation surface. The side wall temperatures are in the range of 1600-1700°C, 
which would cause recession and reduction in the sample diameter. The data from the other two samples are similar; 
stagnation surface temperature exceeded the 3000°C limit in every sample and hence, could not be measured. The 


infra-red (IR) camera also reached its saturation limit at the stagnation surface, therefore, IR data could not be obtained 
for surface temperature. 
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Figure 15: Temperature data from pyrometers and predictions. 


D. Recession measurements 

Pre and post-test shapes were measured by a laser scans in two perpendicular planes (designated as X and Y) 
passing through the geometric axis. Figure 16 shows a representative graphic of an ablated surface (yellow) with 
colored vectors connecting the ablated shape to the initial shape. The scan data from the two planes are collected in 


Figure 17. Side-wall recession reduces the body diameter, which results in tapering of the sides and a decrease in the 
shoulder radius. 
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Figure 16. Laser scan graphic with vectors from initial to Figure 17. Representative pre- and post-test shapes. 
ablated surfaces in two perpendicular Centerline recession is 15.5 mm. 
(X and Y) planes. 


Centerline recession 1s the vertical distance between the ablated and initial surfaces. Table 5 summarizes the results 
from the three PICA tests. Exposure duration and recession varied; however the calculated average recession rate is 
consistent to within 4%. Average recession obtained from CT measurements is 15.62 mm which is consistent with the 
data obtained by laser scans. 

The sample densities also varied, and are all above the nominal PICA value of 0.274 g/cm?. In general, recession 
rate is inversely proportional to density. If the recession data are scaled by initial density, then, as shown in the final 
column of Table 5, the average mass loss rate varies by less than 1%. 

The data From Figure 17 are re-plotted in Figure 18 with the initial shape displaced downward to the ablated 
surface. The latter plot shows that the nose radius is almost unchanged (as expected for an iso-q shape), and the 
diameter decreases by about 7%. 


Table 5. Sample exposure duration and ablation data 


Calculated 
Exposure Centerline Calculated Average Average Mass 
Sample Duration Recession Initial Density Recession Rate Loss Rate 
(+ 0.003 
(+ 0.05 s) (+ 0.25 mm) g/cm3) (mm/s) (g/cm?2-s) 
PICA-1 19.65 15.5 0.278 0.789 0.0219 
PICA-2 20.2 1a 0.28 0.777 0.0218 
PICA-3 20.15 15.3 0.286 0.759 0.0217 
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Figure 18. Comparison of initial and final surface curvatures. 


E. Post-test computational analysis: Test article 

Having gained confidence in the numerical simulation, computations were performed for the test article at the 
heater conditions given in Table 3. Computations for the test article were performed with two different wall boundary 
conditions: (a) wall temperature fixed at 400 K, and (b) wall temperature allowed to float with imposition of radiative 
equilibrium (albeit without conduction through the thickness). Simulations were performed for both the initial and 
final shapes (Figure 17). The predicted surface pressure and heat flux distributions are shown in Figures 19 and 20 for 
the initial and final shapes, respectively. Since the radius of curvature for ablated shape is identical to the initial shape 
of the sample, no significant difference was observed in the heatflux values for ablated shape. However, due to 
recession there is a decrease in static pressure (the distance from nozzle exit plane increases), which reflects in ~5% 
decrease in stagnation pressure as shown in Figure 20. 
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Figure 19. Surface distributions of: (a) pressure, and (b) heat flux from DPLR simulations for the initial shape of test 
article at conditions given in Table 4. Hot-wall computations assume a surface emissivity of 0.85 and heat conduction 
through the thickness of the material is neglected. 
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Figure 20. Surface distributions of: (a) pressure, and (b) heat flux from DPLR simulations for the final (ablated) shape of 
the test article at conditions given in Table 4. The post-test shape is a photographically scaled (0.937-scale) of the original 
shape, and the freestream conditions in the computation corresponding to 15.5 mm movement away from the nozzle exit 
plane. Hot-wall computations assume a surface emissivity of 0.85 and heat condition through the thickness of the material 
is neglected. 


F. Post-test computational analysis: Material thermal response 


FIAT calculations were performed using the nominal hot-wall environments from DPLR and also with +10% 
uncertainty applied to the heat transfer coefficient. Boundary conditions were linearly interpolated in time between 
initial values for the unablated shape (Figure 19) and final values for the fully ablated shape (Figure 20). 


Figure 21 presents the predicted stagnation point recession histories for these three heating levels. The data fall slightly 
below the range of recession predicted by the PICA model. One factor that contributes to the over-prediction is the 
above-nominal virgin density of the MSL-grade PICA used in this work. The data points rise by a few percent if 
corrected for differences in initial density. The standard PICA model appears to be conservative for predicting 
recession at high heat flux. 
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Figure 21. Recession predictions for standard PICA compared with data. 
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VI. Conclusion 


The current arcjet tests on MSL grade PICA material demonstrated the potential of this material at conditions 
beyond the established limits of 1500 W/cm? heat flux and 1.0 atm pressure. The material not only survived but 
performed exceptionally well, with consistent performance across multiple samples, at the estimated heat-flux of 1800 
W/cm and pressure of 130 kPa. While we had problems in obtaining consistent calorimeter data for the heat flux 
values due to melting of slugs, we are confident in the aerothermal computational predictions, as, when same 
aerothermal model was applied to similar heater settings in a past test, it yielded good agreement with 
calorimetry. The predicted stagnation pressure values from aerothermal analysis were in good agreement with 
measured values. The thermal response model (FIAT) predictions for recession show there is slight over prediction 
compared to measured values of recession. More tests at the heat flux ranging from 1500 to 2500 W/cm? as well as 
reliable calorimetry are needed to establish and validate the PICA thermal response model at high heat flux, and 
expand the performance envelope. The present results show a promise to significantly change the design space for 
various missions, where a low density TPS material like PICA could be selected for thermal protection. These missions 
may include sample return missions from various destinations like Mars, comets and asteroids, human space missions 
to Mars as well as science missions to outer planets like Saturn, Uranus and Neptune. 


VII. Future Work 


PICA can be considered as a family of materials rather than a single composition. Stardust, OSIRIS-REx and MSL 
have all flown slightly different variants of PICA and the viable compositional trade space for PICA has not been fully 
explored. The objective for the future test series would be to investigate the performance limits for variants of PICA. 
Based on the success of current arcjet testing for MSL grade PICA, we plan to conduct additional tests with various 
facility settings in order to achieve varying values of heat-flux and pressures. The plan is to achieve heat fluxes in the 
range of 2000 W/cm? - 3500 W/cm”, pressures of 1.5 kPa - 2.0 kPa and shear of 400 Pa - 1000 Pa. All compositions 
evaluated will be compared to baseline MSL PICA. It is anticipated that improvements to key properties of PICA will 
yield a TPS with expanded capabilities and enable future science and human missions. 
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